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Regularity of the Monge— Ampere equation 



^ ■ in Besov's spaces 

c" 

C^ I Alexander V. Kolesnikov and Sergey Yu. Tikhonov 

Abstract. Let fj, = e dx he a probability measure and T = V3> be the 
optimal transportation mapping pushing forward /i onto a log-concave com- 
pactly supported measure u = e^^ dx. In this paper, we introduce a new ap- 
P^ ■ proach to the regularity problem for the corresponding Monge-Ampere equa- 

^ • tion e-^ = detD^* • g-'^'t^*) in the Besov spaces W^^'\ We prove that 

• ^ D^<I> S Wjq'j, provided e~^ belongs to a proper Besov class and W is convex. 

i -~H , In particular, D^^ £ ^foc ^^"^ some p > 1. Our proof does not rely on the 

jrt ' previously known regularity results. 



1. Introduction 



^*7^ , We consider probability measures n — e ^ dx, v — e ^ dx onW^ and the 

' optimal transportation mapping T pushing forward /^ onto v and minimizing the 



Monge-Kantorovich functional 



cn 
en 

O; J \\x-T{x)rdf,. 

It is known (see, e.g., |Vil| . |BoKo2| ) that T has the form T = V<J>, where $ 
is a convex function. If $ is smooth, it satisfies the following change of variables 
formula 

^'. (1) e-^ = e-^(^*'detl?2$. 

C^ ■ 
■ ■ This relation can be considered as a non-linear second order PDE with unknown 

$, the so-called Monge-Ampere equation. 

The regularity problem for the Monge-Ampere equation has a rather long 

history. The pioneering results have been obtained by Alexandrov, Bakelman, 

Pogorelov, Calabi, Yau. The classical theory can be found in |Poj , [Baj , and [GTj . 

See also an interesting survey |Krj on nonlinear PDE's. 
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Despite the long history, the sharpest Holder regularity results of classical type 
have been obtained only in the 90'th by L. Caffarelli |Calj (see also [CCl IGul 
IVilj ). In particular, Caffarelli proved that D^^ is Holder if V and W are Holder 
on bounded sets A and B, where A and B are supports of fi and i/ respectively. In 
addition, B is supposed to be convex. It seems that the latter assumption cannot 
be dropped as demonstrated by famous counterexamples. A nice exposition with 
new simplified proofs and historical overview can be found in [TW] . 

Another result from [Calj establishes sufficient conditions for $ to belong to 
the second Sobolev class W^^'^ with p > 1. More precisely, Caffarelli considered 
solution of the Monge- Ampere equation 

detD'^<P = f 

on a convex set H. with <I>|ao — 0. Assume that fl is normalized: Bi C fl C Bd 
(an arbitrary convex set fl can be normalized using an affine transformation). It is 
shown in [Cal that for every p > there exists £(p) > such that if |/ — 1| < e{p) 
then ||'I*l|vy2,p(B^ ,2) < C{e). Wang jWaj proved that for a fixed e in |/ — 1| < e the 
value of p in the inclusion $ G ^loc cannot be chosen arbitrary large. 

This Sobolev regularity result has been extended and generalized in different 
ways in the recent papers |Sa) . [DPhF] . [DPhFS) . and [Sm'. See also 11] for 
some results on the mean oscillation of D^$. It was shown in DPhFS that every 
$ satisfying det_D^$ = / on a normalized convex set ft with ^\gn — belongs to 

W^'^+^{n'), where n' = (x : $(x) < -i||$||L=o} provided < A < / < A. 

The main purpose of this paper is to develop an alternative approach to the 
regularity problem of the Monge- Ampere equation. We prove that <& belongs to a 
Besov's space under the assumption that e~^ is Besov and W is convex. We give a 
short proof which does not use previously known regularity results. Our estimates 
rely on a generalization of the so-called above-tangent formalism which has been 
widely used in the applications of the optimal transport theory in probability and 
PDE's (see [AGS], |BoKo2| . [ViT] . and [Vi2] V We also apply a result of McCann 
on the change of variables formula and some classical results on equivalence of 
functional norms. 

Some estimates of the type considered in this paper have been previously ob- 
tained in .KolJ in the case of the Sobolev spaces. Applications to the infinite- 
dimensional analysis and convex geometry can be found in [BoKol] and [Ko2j 
respectively. 

Hereafter Br denotes the ball of radius r centered at 0. We use notation D^^ 
for the Hessian matrix of $ and || • || for the standard operator norm. We will 
assume that the measures fi and v satisfy the following assumptions: 

Assumption A: The potential y : M'' — !■ M has the representation 

V^Vo + Vi, 
where |Vb| is globally bounded, Vi admits local Sobolev derivatives, and 

|vyi|eLi(/i). 

Assumption B: The support of j/ is a compact convex set B C Br. 



Definition. Besov's space (Fractional Sobolev's space). The space 
W'P{Q), where Q is a cube in R'*, consists of functions with the finite norm 

iiwik-p(Q) = iimiilp(q) + (y J "J'_^|!+!p dxdyy. 

The space W^-^iR'^) is the completion of C^(R'*) with the norm 

/ f f \u( X) 7/(77 Ir \ ~ 

\\-\W^im-[ll |,_,|.... d.dyy. 

Definition. Log-concave measure. A probability measure i' is called log- 
concave if it satisfies the following inequality for all compact sets A,B: 

v{aA + (1 - a)B) > i^°'{A)i^{B)^-°' 

and any < a < 1. If i^ has a density v = e~^ dx, then W must be convex (we 
assume that W — +00 outside of supp(i^)). This is a classical result of C. Borell 



o|. 

Remark 1 . We note that the second derivative of the convex function $ can 
be understood in different ways. In the generalized (weak) sense this is a measure 
with absolutely continuous part D^^ dx and singular part D^^. Throughout the 
paper we use the following agreement: the statement "_D^$ belongs to a certain 
Besov or Sobolev class" means that the measure D^^ has no singular component 
and the corresponding Sobolev derivative Z3^$ = D^^ belongs to this class. 

Theorem 2. Let Assumptions A and B be fulfilled and, moreover, 
(1) there exist p > I and < 7 < 1 such that for every r > 



Br \yf+^^'' 



dy < 00, 



where SyV = V{x + y) + V{x — y) — 2V{x) and {5yV)+ is the non-negative 
part of 6yV ; 

(2) e^^ is locally bounded from below; 

(3) v is a compactly supported log-concave measure. 

Then 

\\D^Hw^-HQ) < °° 
for every cube Q C M.'^ and < £ < ^ . 

Taking p = +cx) we obtain the following result. 

Corollary 3. Let Assumptions A, B and conditions (2)- (3) he fulfilled. As- 
sume, in addition, that {5yV)^. < w(|y|) for some uj : M^ — )■ M^ and 

u}{s) 



I 

Jo 



. , ds < 00 

1+7 
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for some 7 > and every r > 0. Then |jZ? $||vi/e.i(q) < cxd with < £ < ^ 

In particular, applying the fractional Sobolev embedding theorem (see, e.g., 
[Adl Ch. V]) we get 

Corollary 4. Under assumptions of Theorem\^ for every e > 0, 



\D'n^Li--^'- 



Remark 5. Note that in corollary |3] we do not assume that V is bounded from 
below. 



2. Auxiliary results 

Below we will use the following function space (see, e.g., [St]). 
Let A^'', < a < 2, be the space of functions with the finite norm 

{\\f{x + t) + f{x~t)~2f{x)\\py 



l/ll 



\t\d+aq 



dt 



where || • \\p is the LP-norm with respect to the Lebesgue measure. We will apply 
the following equivalence result from [St) Ch.5, Sec. 5, Propos. 8']. 

Lemma 6. For a > 1 the norm ||/||a''''' is equivalent to 

(||V/(x + t)-V/(a:)y^ ,, 

\t\d+{a-l)q 

In particular, if a > 1 and ||/||Ag''' < oo, then f admits the Sobolev derivatives. 

Let us recall that every convex function $ on R'' admits the following two types 
of second derivatives. 

Definition. We say that the measure fj,ev is the distributional derivative of 
a convex function $ along unit vectors e, v if the following integration by parts 
formula holds for any test function ^: 

£, dflev ^ - I deS, dv^ dx. 

We set 

dev^ ■= fJ.ev 

Definition. The absolutely continuous part {dev^)a of dev^ is called the 
second Alexandrov derivative of $ along e,v. Clearly, 

See* > (9ee*)a > 

in the sense of measures. 

Let us denote by D^^^ the matrix consisting on these absolutely continuous 
parts. We will apply the following result of R. McCann from [McCannJ. 

Theorem 7. For ^-almost allx the following change of variables formula holds 

e-^(-)=detZ?2<i>(^).e-M'(v*(-)). 

We will need the following lemma from |Kolj . In fact, this is a generaliza- 
tion of the well-known above-tangent lemma which has numerous applications in 
probability and gradient flows of measures with respect to the Kantorovich metric 
(see [ViT] . [Vi2] . [AGS] , |BoKo2j ). The proof follows directly from the change of 
variables and integration by parts. 



Lemma 8. Assume that W is twice continuously differentiable and 

D^W >K -Id 
for some K ^R, p > 0. Then 

SyV{dy<S>)Pdn 
- f / 1^*^^ + ^^ ~ V$(x)p(^y$)P d^i+^ f |V$(x - y) - V$(x)p(5j,$)P d^i 

where SyV = V{x + y) + V{x - y) - 2V{x). 

Corollary 9. It follows easily from Lemma\^that inequality 

J 5yV{5y^fd^l >pf (V5y^, ( i^^ $ ) " 1 V J^ $ ^ ( J^ $ )P- 1 dfl 

holds for any log-concave measure v. In particular, this holds for the restriction of 
Lebesgue measure tt^-^U '''^ ^ convex subset A. In this case W is a constant on 

A and W{x) = +00 if x ^ A. 

Remark 10. Let us assume that V is twice differentiable and y = te for some 
unit vector e. Dividing by t^P+2 and passing to the hmit we obtain 

(2) l^ee'&LrfM > K j \\D^^ ■ ef^l, dp. 

+ p /"((I?2$)-lV$ee,V$ee)$L"'rfM• 
N0W it is easy to get some of the results of [Kolj from ([2]). In particular, 
applying integration by parts for the left-hand side and Holder inequalities, one 
can easily obtain that 

^ll*eellL.(p)<^llK'llL.(^)- 

It is worth mentioning that ^ gives, in fact, an a priori estimate for derivatives 
of $ up to the third order (due to the term p/((D2<I>)~iV$ee, V<i>ee)*ecr^ dp)- 

We denote by Ao$ the absolutely continuous part of the distributional Lapla- 
cian of $. 

Lemma 11. Under Assumptions A and B, there exists C such that 

I Aa$(x + y)dp<Y, f e'""'"^ d[d,^,Mx + v)] < C 

uniformly in y ^ M."^. 

Proof. The inequality / Aa$(a; + y) dn< J^t^i I e^^^^^ d[dx,x,^{x + y)] is 
clear in view of the fact that the singular part of (9a;ia;i<I' is nonnegative. Moreover, 
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we get 



= ci /"(V$(x + y),VFi(x))e-^i(^) dx < ciR /" |VFi(x)|e-^i(^) dx 
<C2 / |VV^i(a;)| d^l. 



D 

Finally, we will use the fractional Sobolev embedding theorem (see |Ad|, Ch. 
V]). We formulate it in the following form given in |MSh| (see also [BBM] . |KL) ). 

Theorem 12. Letp > 1, < s < 1 and sp < d. Then for every u e Wq'^{R'^) 
one has 

ll'"llL9(K<i) - c(d,p) .^ _ \p-i ll"llM/g=-P(M<i)' 

where q — dp/{d — sp). 



3. Proof of Theorem [2] 

Let us apply Lemma [8] with p = 1. We have 

/ SyV-Sy<^dH> I {V5y^,{Dl^)-^V5y^) d^l. 

Taking into account Lemma [TTl we obtain J^^ ll^a*^!! d/i < oo. Then Cauchy 
inequality yields 

\\Dl<i>\\dfl- I SyV-Sy^dH>(f \\/Sy'i>\d^y . 

By the Holder inequality, for every p, <7>1, — h- = l, we have 

byV ■ Sy^ df,< f {SyV)+ ■ Sy^ dfi < 1 1 (5, F) + | ^P (^) • 1 1 <5, $ | ^ , (,,) 

JR'' 

Let us now estimate ||(5y$||ig(^). Note that |V$| < R, and hence 6y^ < 2R\y\. 

Let us also mention that t — > ^(x + ty) is a one-dimensional convex function for 
a fixed a;. Therefore, ruy = dff^{x + ty) is a nonnegative measure on M}. Moreover, 

Sy^ = / (V$(a; + sy) - V$(x - sy),y) ds 
Jo 

= / dmy ds. 

Therefore, we have 

W^ynU^^ ^ J {Sy^r dfi < {2R\y\y~' J dy<^ dp 

= (2R\y\Y'^ j f j e-'^^'-'U[dl^{x + ty)\dtds. 



/O J-s. 
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It follows from Lemma [TT] that 



WSynU^^SUiy 



<C|w|l+«. 



Hence, 



(4) 



JR'' ' Jm'' 

Now we divide this inequality by |2/|''^^^''' and integrate it over a bounded subset 
Q cM.'^. By condition (1), we get 

(3) / rjp^(y"jV$(x + y) + V$(a;-y)-2V$(a;)|d/i)'d2;<^. 

Let us take a smooth compactly supported function ^ > 0. We will show that 

I J+.+j f \^{x+y)S7^x+y)+^ix~y)\7^x-y)-2^ix)V<^ix)\ dxY dy< oc. 

nd \y\ ^-^ ' ^jRd / 

Let us split this integral in two parts: 

/ ■ • ■ dx ^ ■ ■ ■ dx + ■ ■ ■ dx = Ii + I2 ■ 

JRrf JBi JBf 

To estimate the second part, we note that 

/ \^{x + y)\7^{x + y) + ^{x-y)y^{x-y)-2^{x)V<^{x)\ dx<4R \^{x)\ dx. 

Thus, I2 < 00. 

Let us estimate /i. It follows from estimate ([3|) and condition (3) of the theorem 
that 

/ I H+2+A l |V$(x + y) + V$(x~2;)-2V$(a:)|C(a;)dx)'dy<oo. 

JBi \y\ VH<i ^ 

Therefore, it is enough to show that 

X 2 

+ V^{x-y){£,{x-y)-£,{x))\ dx\ dy < 00. 
Since |y| < 1 and ^ is compactly supported, there exists Rq > such that 
1 



/:^< 



Bno \y 



d+2+1 



\V<i>ix + y){i{x + y)~C{x)) 



+ V<Pix~y){^{x~y)~^{x))\dx^ dy < 
Using smoothness conditions on ^, we obtain 

V$(a; + 2/)(e(x + y) - ^{x)) + V$(x - y){C{x ~ y) ^ ^{x)) 

= (V$(x + y)-V$(a;-y))(e(a; + 2/)-C(a^)) 

+ V$(a:; - y) {^{x + y) + ^{x - y) - 2^{x)) 

< c(i?o|y|' + \y\\V^x + y)- V$(x - 2/)|^ 



Thus, it is sufficient to show that 

(5) f Y^il \V<^{x + y)~y^[x~y)\dx) dy < ^. 

To prove ([S]), we use the representation 

/ |V$(x + y)- V$(x-y)| da;= / / S\d\ds^eX^ ^ ^y) {s) ■ y^\ dx 

J Br J Br J -I ._i -' 



i=l 



The latter is bounded by C\y\ J„ A$. This immediately implies that I^ < oo and 
therefore @ is proved. 

This means that |^ • V<i>|^i,2 < cx) for every smooth compactly supported ^. 

Using smoothness conditions on ^ and boundedness of V$, we get from Lemma [6] 
that D^^ has no singular parts and 



(6) 



L W+^^L l^''^(^ + 2^)-^''^(^)l^^)''^2/<°° 



for any B^. 

Finally, applying the Cauchy-Schwarz inequality, we obtain for every 5 > Q 

Changing variables implies 



dx dz < 00 

QJQ 



\z-x\'^+'' 



for every < e < ^ and bounded Q. The proof is now complete. D 



4. Remarks on improved integrability 

d 

By applying Theorem[2]we get a better (local) integrability of ||D'^$|| (L''" 2+^ 
instead of i^); see CoroUarylH This can be used to improve the estimates obtained 
in Theorem [2] 

We assume for simplicity that we transport measures with periodical densities 
by a periodical optimal mapping 

T[x) = .T + V<y5(a;), 

where Lp is periodical. Equivalently, one can consider optimal transportation of 
probability measures on the flat torus T. Then one can repeat the above arguments 
and obtain the same estimates which become global. 

In general, it can be shown that the assumption ||-D'^(y5|| £ L'^iT) implies 

D'^u){x + y) - D'^Lp{x)\'^ , , 

- ax dy < 00 



JTJJ Ij/I'^+tttCt+V)-'^ 



for every e and ||-D^<y3|| G U' (T) with any r' satisfying 



r < 



2dr 

r+1 



- 7fT(7 + ^) 

Starting with rg = 1 and iterating this process one can obtain a sequence r„ such 
that ||I?^(/3|| e U"^^ for every r„ and e > 

2dr„ 
r^i + l ^ ' q ' 

One has ||I?^(^|| G U'~^ ^ where r = Hni„ r„ solves the equation 

x^ + a;(q(7 - d) - 1) + qd = 
with r > 1. 
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